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Abstract 



We review the production of high-mass lepton pairs in fixed-target experiments, 
including both Drell-Yan (DY) and heavy quarkonium production (J/ip, Y>', T(IS'), 
T(2S), and T(3S)). In recent years DY data has become standard input to the de- 
termination of parton density distributions. DY data has recently yielded the first 
measurement of the x-dependence of the d/u asymmetry of the proton. Similar to 
the observations in deeply inelastic scattering, precision measurements of the nuclear 
dependence of the proton-induced DY process exhibit shadowing at small target mo- 
mentum fraction, xi- There is, however, no evidence of enhanced DY production from 
nuclear targets. Mean transverse momenta of DY pairs is observed to increase with 
target mass. These data, analyzed within a new theoretical framework, provides an 
estimation of the energy loss of fast quarks in nuclear matter. In contrast to the DY 
process, there are large nuclear effects in the production of all quarkonia. These effects 
show strong dependence on p t and xf, but do not scale with X2- Statistically less 
precise measurements of D meson (open-charm) production at small xf are consistent 
with no nuclear dependence. 
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1. INTRODUCTION 

The experimental detection of high-mas^] lepton pairs produced in hadronic reactions 
has a long and rich history. The famous quarkonium states that revealed the existence of 
the charm and beauty quarks in the 1970s were discovered through their dilepton decay 
branches. They are superimposed on a continuum, which was anticipated theoretically in 
1970 |L|, and is now known as the Drell-Yan (DY) process (Figure 0). The DY process, 
electromagnetic quark-antiquark annihilation, is shown diagramatically in Figure along 
with its close relative, deeply inelastic lepton scattering (DIS). By 1980, DY production was 
already a source of information about antiquark structure of the nucleon. Additionally, DY 
production with beams of pions and kaons yielded the structure functions of these unstable 
particles for the first time. Also notable in the history of the DY process were the discoveries 
of the W ± and Z° particles in 1983, produced by a generalized (vector boson exchange) 
quark-antiquark annihilation mechanism. 

In the mid 1980s some nagging theoretical issues were resolved, with the result that by 
the 1990s the DY process had joined DIS as a source of quantitative information on the 
quark structure of hadrons. Because the DY process at leading order explicitly involves 
one antiquark, it provides a sensitive measurement of antiquark structure even in kinematic 
regions where quarks dominate in number. 

The peaks in Figure [I] are interlopers in an otherwise well understood spectrum. They are 
seen along with the DY pairs because the 3 Si states of heavy quarkonium have substantial 
decay branches to dileptons. However, the production of heavy quarkonium is a strong 
interaction process. At beam energies well above threshold, gluon fusion is the dominant 
production mechanism. Unlike the well-characterized DY process, the theoretical description 
of quarkonium production at the low transverse momenta (p t ) characteristic of fixed-target 
experiments requires a number of model-dependent assumptions . 

Although the production mechanisms of quarkonia and the DY continuum differ greatly, 
these two types of processes are often reviewed together ||, [|, H, (Tj because the data come from 
the same experiments. We continue that tradition here. The main emphasis of this review is 
the application of dilepton production as a tool to investigate particular physics topics. New 
experimental work has been carried out in recent years by few but prolific collaborations 
working in the fixed-target programs at the CERN SPS accelerator and at Fermilab. This 
review concentrates on new experimental results published after the comprehensive review 
of Freudenreich |7j . 

Lepton-pair production forms a major thrust of the collider program at Fermilab. We 
review those results only as they pertain directly to the physics addressed by the fixed-target 
experiments. 

Section 2 presents an overview of the DY process, highlighting significant developments 
in theory and experiment. Sections 3.1 through 3.4 present recent and future applications of 
the DY process to the study of the quark structure of hadrons. In the remainder of Section 
3, we discuss several physics topics related to the nuclear dependence of the DY process. 

1 In the context of this review, high mass is M;+;- > Mj/^ 
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Section 4 reviews the status of quarkonium and open heavy-flavor production in nuclei and 
their relation to quarkonium suppression as a signal for the formation of the quark-gluon 
plasma (QGP) in high-energy collisions of heavy ions. Finally, in Section 5, we look at some 
of the near-future opportunities in dilepton physics. 

2. THE DRELL-YAN PROCESS 

2.1 Drell-Yan Process in the Parton Model 

The parton model expression for the DY process conveys the essential simplicity of the 
reaction. 

d^~ F = Imh ( Xl I x 2 ) ? + (i) 

Here q a (x) are the quark or antiquark structure functions of the two colliding hadrons eval- 
uated at momentum fractions x\ and x%- The sum is over quark flavors. In addition, one 
has the kinematic relations, 

r = X\X 2 = M 2 /s, 

x F = Xx- x 2 , (2) 

where M is the invariant mass of the lepton pair and s is the square of the center-of-mass 
energy. The cross section is proportional to a 2 , indicating its electromagnetic character. The 
parton-model scaling properties of the DY process are illustrated in Figure |3], which shows 
data from experiments at 400 GeV/c M and at 800 GeV/c (|)], TO] ; PL McGaughey et al, 
unpublished data). 

In the parton model, the angular distribution of dileptons is characteristic of the decay 
of a transversely polarized virtual photon, 

^ = a o (l + W0), (3) 

where 9 is the polar angle of the lepton in the virtual photon rest frame and A = 1. Early 
experimental data from both pion and proton beams |3J] were consistent with this form but 
had large statistical errors. 

Recently, E772 has performed a high-statistics study of the angular distribution for 



DY events |T2| with masses above the T family of resonances. About 50,000 events were 
recorded from 800 GeV/c p + Cu collisions, using a copper beam dump as the target. 
Figure |] shows the acceptance-corrected angular distribution, integrated over the kine- 



matic variables. Analyzed in the Collins-Soper reference frame |0|, the data yield A = 
0.96 ± 0.04 ± 0.06(systematic). 

2.2 QCD and the Drell-Yan Process 

Logarithmic scaling violation of the DY cross section is expected in leading order as 
for DIS. In the DY process, however, the virtual photon is timelike, M 2 = Q 2 > 0. Di- 
rect experimental confirmation has been difficult for the DY process (see the discussion in 
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Freudenreich's review 0) for two reasons. First, there is the complication that the cross 
section involves the convolution of two structure functions. Second, and more importantly, 
the DY process is experimentally well established as the major contribution to the dilepton 
spectrum only for M > 4 GeV/c 2 . This is already a large value of Q 2 by DIS standards. 
Increasing M by a factor of five or so to provide enough lever arm to see the logarithmic 
change in q(x, M 2 ) involves measurement of an exceedingly small cross section. Now that 
DY data has become an integral feature of global parton-structure fitting programs, there is 
no longer much concern about testing its logarithmic scaling features. An excellent summary 
of the use of DY data in structure function phenomenology appears in a recent review fll4 



The importance of next-to-leading order (NLO; terms proportional to a s ) QCD cor- 
rections to the DY process has been well known since the mid-1980s || |15| . Historically, 
two experimental features demanded theoretical improvement: first, the experimental cross 
section was about a factor of two larger than the parton-model value, and second, the distri- 
bution of dilepton transverse momenta extended to much larger values than are characteristic 
of the convolution of intrinsic parton momenta. Corrections to the DY mechanism at NLO 
largely fix these two problems. As shown in Figure |5], this comes at the expense of involving 
diagrams that seem to eliminate the simple parton interpretation of qq annihilation. Nev- 
ertheless, for most interpretations of the DY process, it is reasonably accurate to think in 
terms of parton-level qq annihilation with a modest increase in cross section. 

QCD provides a procedure for calculating higher-order corrections to the DY process 
as well In practice, NLO corrections are normally employed in QCD phenomenology. 
Figure |3| shows the typical agreement found when NLO calculations are compared — without 
normalization — to data. 



Including higher-order QCD corrections to the DY process [17], [18| results in the more 
complicated form of the angular distribution, 

oc 1 + \cos 2 9 + nsin29cos(fi H — sin 2 9cos2(f), (4) 
ail 2 

where <ft is the azimuthal angle and A, fi, and v are angle-independent parameters. NLO 
calculations predict small deviations from 1 + cos 2 9 (< 5%) for p t below 3 GeV/c. The 
relevant scaling parameter for the magnitude of these deviations is Pt/Q, implying that NLO 
corrections become important when p t — Q. A relation, 1 — A — 2u = 0, developed by Lam & 
Tung 1 201 , is analogous to the Callan-Gross relation in DIS. Measurements with pion beams 



at CERN jnj] and at Fermilab JEJ have shown that the Lam- Tung relation is clearly violated 
at large p t . 

Pion-induced DY experiments have unexpectedly shown that transverse photon polar- 
ization changes to longitudinal (A ~ —1) at large xp [^1], |22|, |24J. The xp dependence 



of A is qualitatively consistent with a higher-twist model originally proposed by Berger & 
Brodsky pR p6| . However, the quantitative agreement is poor. The model's basis can be 
described as follows. As xp of the muon pair approaches unity, the Bjorken-x (momentum 
fraction) of the annihilating projectile parton must also be near unity. Thus, the whole pion 
contributes to the DY process. This can be treated with perturbation theory, with the result 
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that the transverse polarization of the virtual photon becomes longitudinal. The angular 
distribution at large xp becomes 



^ oc (l-x) 2 (l + \cos 2 6) + asin 2 6, (5) 

where a is oc p 2 /Q 2 - 

Eskola et al |27| have shown that an improved treatment of the effects of nonasymptotic 
kinematics greatly improves quantitative agreement with the A values from the pion data. 
Brandenburg et al |28| have extended the higher twist model to specifically include pion 



bound-state effects. They predict values for A, fi and v that are in good agreement with the 
pion data at large xp. Unfortunately, the results are quite sensitive to the choice of the pion 
Fock state wave functions, which are not well constrained by experimental data. 

At this time, no comprehensive theory covering the full kinematic range of the experi- 
ments is available, nor is the observed violation of the Lam- Tung rule understood. 

2.3 Factorization and Nuclear Size Effects 

General discussions of factorization in the DY process and other QCD processes may be 



found in reviews by Collins & Soper |29] and Sterman et al 

A crucial requirement for the application of the DY process for the study of the quark 
structure of nuclei and nucleons (where most targets are actually nuclei) is that the factorized 
form of the cross section expressed by Equation [I] (or its QCD-corrected equivalent) remains 
valid for nuclei. Specifically, one wants assurance that q a {x, M 2 ) is the same for a proton that 
has traversed the diameter of an atomic nucleus as it is for a free proton. Fortunately, this 
problem was analyzed extensively in the 1980s, particularly for the DY process |31~| , |32"| , p3| . 
Within limits given approximately by 

Q 2 > x 2 m N L A (l 2 L ), (6) 

for hadron-nucleus collisions, initial-state interaction effects can be ignored for observables 
integrated over dilepton transverse momenta. In Equation [| x 2 is the parton momentum 
fraction for a nucleus of mass A, m N is the nucleon mass, and Zj_ ~ 0.3 GeV/c is the average 
momentum exchange in a quark-nucleon collision. With a path length La = ZRqA 1 ^, 
Equation |6| becomes Q 2 > X2A 1 / 3 GeV 2 /c 2 , a condition readily met for all nuclei. New nuclear 
dependence measurements, discussed in Sections 3.5-3.7, provide quantitative evidence of the 
absence of substantial initial-state interaction effects. 

2.4 Experimental Detection of Continuum Lepton Pairs 

The DY process was anticipated theoretically in 1970 (]]]. However quantitative ex- 
periments had to wait until the late 1970s for the development of techniques permitting the 
measurement of picobarn cross sections in the presence of background processes many orders 
of magnitude larger. The most successful high-luminosity DY experiments have exploited 
inclusive dimuon detection in magnetic spectrometers whose apertures are filled with hadron 
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absorbers. Hadron absorption near the target also reduces the backgrounds of muons from 
pion and kaon decays. Freudenreich's 1990 review [[?[] of lepton-pair production includes 
detailed descriptions of spectrometers built and operated during the 1970s and early 1980s. 
During the late 1980s and 1990s, two dimuon spectrometers remained operational at Fermi- 
lab and two at the CERN SPS. This review concentrates on data from these instruments. 
(The HELIOS collaboration at the SPS focused on the mass region below the J/ip; hence, 
their work is not covered by this review.) 

At Fermilab, a series of three experiments — E772, E789, and E866 — used differing 
configurations of the large spectrometer located on a direct 800-GeV/c proton beam line in 
the Meson-East experimental area. The spectrometer (Figure ^) was originally constructed 
for E605, which studied dimuon production for M„+„- > 7 GeV/c 2 . The first two dipole 
magnets focus high-mass pairs into the spectrometer acceptance, thus avoiding the beam 
dump contained in the second magnet. The dump absorbs not only the beam but also 
the enormous flux of low-pt secondaries from the target. The third large analyzing magnet 
provides a p t kick of about 1 GeV/c for momentum analysis by the tracking system. The 
spectrometer is capable of operating at luminosities in the range 10 35 cm _2 sec _1 with 10% 
interaction length targets and detection acceptances of ~ 1 — 2% for the DY process. E772, 



E789, and E866 also exploited the copper beam dump as a second target [12, |3J. The beam 



dump "target" has the advantage of a more open acceptance, permitting the measurement 
of angular distributions of muon pairs over wide range of decay angles. 

The E705 collaboration at Fermilab operated an open-aperture (no hadron absorber) 
spectrometer which contributed important data to the study of the J/ip and ip' resonances 
and their decays f35 |. 

The muon-pair spectrometer operated most recently at CERN was built originally for 
experiment NA10 | 36fl . It consists of a large air-core toroid magnet with tracking chambers 
following the hadron absorber. In recent years, with substantial changes in its configuration, 
it has made unique contributions to the study of J/ip production in relativistic heavy-ion 
collisions at the SPS P7| , as well as to the determination of the d/u asymmetry of the 



proton [38 



2.5 Backgrounds in Inclusive Dimuon Experiments 

The E866 spectrum (Figure |1]) illustrates the typical experimentally "safe" regions for DY 
muon pairs in fixed-target experiments with beams in the few hundred GeV/c range. Below 
the J/ip , a number of potential backgrounds, including the semileptonic decay of charmed 
hadrons, make the interpretation of the continuum more complex. Additionally, depending 
on details of the target and absorber configuration, pion and kaon decays can affect the 
low-mass spectrum. The data can be corrected for the latter contribution via subtraction 
of properly normalized like-sign muon pairs. But charm pair decay, where this technique 
does not work, presents a more serious problem and ultimately limits the kinematic region 
that can be exploited. Ideally, the charm-pair decay contribution could be separated by 
use of a microvertex detector. To date, however, such a detector has not been combined 
with closed-aperture dimuon detection. The most systematic study of the contribution of 
charm production to the dimuon spectrum has been performed recently by the NA50 col- 
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laboration |37| in connection with measurements of J / ip suppression in heavy- ion collisions. 
In NA50, the charm-pair decay continuum is evaluated with the help of simulations of all 
contributing processes in the mass region just below the J/ip. 

3. DRELL-YAN PROCESS AND PARTON DISTRIBUTIONS 

3.1 d/u(x) Asymmetry of the Proton 

From neutrino-induced DIS experiments, it is known that the strange-quark sea in the 



nucleon is roughly a factor of two less than the up or down quark sea [p9|| . The lack of SU(3) 
flavor symmetry in the nucleon sea has been attributed to the much heavier mass of the 
strange quark. Until recently, it had been assumed that the distributions of u and d quarks 
were identical. Although the equality of u and d in the proton is not required by any known 
symmetry, it is a plausible assumption for sea quarks generated by gluon splitting. Because 
the masses of the up and down quarks are small compared to the confinement scale, nearly 
equal numbers of up and down sea quarks should result. 

The assumption of u(x) = d(x) can be tested by measurements of the Gottfried inte- 
gral [HO], defined as 



i r 



where F$ and K 



F%(x,Q 2 )-F%{x,Q 2 ) fx dx 



1 2 

3 + 3 



1 r 



u p \x i 



- d p (x) 



dx. 



(7) 



2 auLi ± 2 are the proton and neutron structure functions measured in DIS ex- 
periments. The second step in Equation [7| follows from the assumption of nucleon charge 
symmetry (Section 3.3). Under the assumption of a symmetric sea, u = d, the Gottfried 
Sum Rule (GSR) 



HJ, Iq — 1/3, is obtained. Several early attempts |4l], |42|, |43fl to test the 
GSR were inconclusive because the measurements did not reach small enough x, where large 
contributions to the integral are expected. Nevertheless, Field & Feynman interpreted 
the early SLAC data [gTJ, which gave J G (0.02 -> 0.8) = 0.20 ± 0.04, as an indication that 
the GSR is violated with d(x) > u(x). They suggested that the Pauli-blocking effect sup- 
pressed g — > uu relative to g — > dd, since protons contain two w-valence quarks and only one 
o?-valence quark. 

The most accurate test of the GSR was reported in 1991 by the New Muon Collaboration 
(NMC) |[45|| , which measured Ff and F% over the region 0.004 < x < 0.8. They determined 
the Gottfried integral to be 0.235 ±0.026, significantly below 1/3. This surprising result has 
generated much interest. 

Although the violation of the GSR can be explained by assuming unusual behavior of 
the parton distributions at very small x [JPJ], a more natural explanation is to abandon the 
assumption u = d. Specifically, the NMC result implies 



d(x) 



u[x 



dx = 0.148 ±0.039. 



We emphasize that only the integral of d — u was deduced from the DIS measurements. The 
x dependence of d — u remained unspecified. 

The proton-induced DY process provides an independent means to probe the flavor asym- 
metry of the nucleon sea ||47||. An important advantage of the DY process is that the x 
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dependence of d/u can be determined. It is interesting to note that, as early as 1981, Fermi- 



lab E288 [48| reported evidence for a d/u asymmetry, based on a measurement of the p + d 
DY cross section. However, this interpretation depended sensitively on assumptions about 
the shape of the valence quark distributions and was not conclusive. Later, the Fermilab 
E772 collaboration |^9| compared the DY yields from isoscalar targets with that from a 
neutron-rich (tungsten) target, and constraints on the nonequality of u and d in the range 



0.04 < x < 0.27 were set. More recently, the CERN experiment NA51 [|38| carried out a 
comparison of the DY muon pair yield from hydrogen and deuterium using a 450 GeV/c 
proton beam. They found that u/d = 0.51 ± 0.04 ± 0.05 at (x) = 0.18, a surprisingly large 
difference between the u and d. 

A DY experiment (E866), aiming at higher statistical accuracy and wider kinematic 
coverage than NA51, was recently completed || at Fermilab. This experiment also measured 
the DY muon pairs from 800-GeV/c protons interacting with liquid deuterium and hydrogen 
targets. The acceptance of the spectrometer was largest for xf = x\ — x 2 > 0. In this 
kinematic regime the DY cross section is dominated by the annihilation of a beam quark 
with a target antiquark. The DY cross section ratio at large xp is approximately given as 

(ypvip + d) d(x 2 ) \ 

2a DY (p + p) ~ 2 ^ u(x 2 )J ' 1 ) 

The ratio is unity when d = u. Figure [j] shows that the E866 measurement of this ratio 
clearly exceeds unity for an appreciable range in x 2 . 

Using an iterative procedure ||, values for d/u were extracted by the E866 collaboration 
at Q 2 = 54 GeV 2 /c 2 . These are shown in Figure |8] along with the NA51 measurement. For 
x < 0.15, d/u increases linearly with x and is in good agreement with the CTEQ4M |TTJ and 
MRS(R2) |5(J parameterizations. However, a distinct feature of the data, not seen in either 
parameterization, is the rapid decrease toward unity of d/u beyond x = 0.2. 

The d/u ratio, along with the CTEQ4M values for d + u, was used to obtain d — u over 
the region 0.02 < x < 0.345 (Figure H). Being a flavor nonsinglet quantity, d(x) — u(x) 
is decoupled from gluon splitting. From the results shown in Figure || one can obtain an 
independent determination |5l| of the integral of Equation |8|. E866 finds 0.100±0.007±0.017, 



consistent with, but roughly 2/3 of the value deduced by NMC. 

Very recently, the HERMES collaboration reported a semi-inclusive DIS measurement of 
charged pions from hydrogen and deuterium targets |52fl . Based on the differences between 
charged-pion yields from the two targets, the ratio (d — u)/(u — d) is determined in the 
kinematic range, 0.02 < x < 0.3 and 1 GeV 2 /c 2 < Q 2 < 10 GeV 2 /c 2 . The HERMES results 
for d — u, shown in Figure [|, are consistent with the E866 results obtained at significantly 
higher Q 2 . 

The E866 data clearly affect the current parameterization of the nucleon sea. The most 
recent structure functions of Martin et al |53| (MRST) included the E866 data in their global 
fit and are very different from the previous MRS(R2) parameterization (Figure^). What is 
not so obvious is that the E866 data also affect the parameterization of the valence-quark 
distributions. Figure |TU| shows the NMC data for F^—F^ at Q 2 = 4 GeV 2 /c 2 , along with the 
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fits of MRS(R2) and MRST. It is instructive to decompose F$' 
from valence and sea quarks: 



x 



x 



F?(x) 



1 r 

-x [u x 



(x) — d v (x)] + —x u(x) — d(x) 



(x) into contributions 



(10) 



As shown in Figure [10], the E866 data provide a direct determination of the sea-quark 
contribution to F^ — F^. In order to preserve the fit to F% — F£, the MRST parameterization 
of u v — d v is significantly different from that of MRS(R2). Neither set reproduces Ff — Fg 
in the range x = 0.2 — 0.3. 



We now turn to the origin of the d/u asymmetry [51, 54]. As early as 1983, Thomas |55| 
pointed out that the virtual pions that dress the proton will lead to an enhancement of d 
relative to u via the (nonperturbative) "Sullivan process." Sullivan |5£]] previously showed 
that in DIS virtual mesons scale in the Bjorken limit and contribute to the nucleon structure 
function. Following the publication of the NMC result, many papers |)7], [ 



59, 60, 61, m 



63. 64, 65] treated virtual mesons as the origin of the d/u asymmetry (see [54] for a recent 
review). Here the n + (du) cloud, dominant in the process p — > 7T + n, leads to an excess of d 
sea. 

A different approach for including the effects of virtual mesons has been presented by 
Eichten et al |63| and further investigated by other authors |65|. In chiral perturbation 
theory, the relevant degrees of freedom are constituent quarks, gluons, and Goldstone bosons. 
In this model, a portion of the sea comes from the couplings of Goldstone bosons to the 
constituent quarks, such as u — > dn + and d — ► ui\~ . The excess of d over u is then simply 
due to the additional valence u quark in the proton. 

The x dependences of d — u and d/u obtained by E866 provide important constraints 
for theoretical models. Figure |9| compares d(x) — u(x) from E866 with a virtual-pion model 
calculation, following the procedure detailed by Kumano |58|. A dipole form, with A = 1.0 
GeV for the tcNN form factor and A = 0.8 GeV for the 7riVA form factor, was used. A 
is the cutoff parameter for the pion form factor. Figure |] (dotted curve) also shows the 
predicted d — u from the chiral model. We follow the formulation of Szczurek et al |55| 
to calculate d(x) — u(x) at Q 2 = 0.25 GeV 2 /c 2 and then evolve the results to Q 2 = 54 
GeV 2 /c 2 . The chiral model places more strength at low x than does the virtual-pion model. 
This difference reflects the fact that the pions are softer in the chiral model, since they are 
coupled to constituent quarks, that carry only a fraction of the nucleon momentum. The x 
dependence of the E866 data favors the virtual-pion model over the chiral model, suggesting 
that correlations between the chiral constituents should be taken into account. 

The chiral and the meson-cloud models both predict that the u and d quarks will carry 
negligible amount of the proton's spin |ET], [63|, B6| , a prediction that remains to be tested. 



Another interesting consequence of the meson-cloud model is that the s and s distributions 
in the proton could have very different shapes, even though the net amount of strangeness 
in the proton vanishes |)7], |68| . By comparing the v and v induced charm production, the 
CCFR collaboration found that the s and s distributions are very similar [p9] |. Dimuon 
production experiments using K ± beams might provide an independent determination of 
the s/s ratio of the proton. 
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3.2 Polarized Drell-Yan and Production 



The spin structure of the nucleon has been extensively investigated by polarized DIS [70 



7TJ during the past 10 years. Much of the excitement in this field has come from the finding 
that quarks contribute a surprisingly small fraction to the nucleon's total helicity. As a 
result, considerable interest has been centered on QCD processes in hadronic collisions, 
where polarized gluon and antiquark effects are directly manifested. 

Polarized DY production and W ± production have great potential for providing quali- 
tatively new information about antiquark polarization.^ Because polarized proton collisions 
will become feasible at the RHIC facility starting in 2001, it is timely to review the type of 
information that such studies will provide. 

The longitudinal spin asymmetry (requiring polarization of both protons) in the DY 
process is, in leading order, given by 

[Aq a (x 1 )Aq a (x 2 ) + Aq a (x 1 )Aq a (x 2 )} , , 

' LL " l " 2 Y,a e2 a\ ( la{x l )q a (x2) + qa{x l )q a {x 2 )\ 

with Aq a = g+ — q~. The superscripts refer to parton spin projections parallel (+) or 
antiparallel (— ) to the proton's spin projection.^ One can simplify this equation by choosing 
the kinematic region xp > 0.2, where a single term, u(xi)u(x 2 ), dominates the denominator 
of Equation |TT] (See Section 3.1). The equation thus becomes 

/\Ti 

A%(x u x 2 ) » < S (^i) x -fe), (12) 

u 

where the asymmetry measured in DIS on polarized protons is A p L l L IS {x) = g\(x)/F\(x). 
Thus, the helicity asymmetry in polarized p-p collisions is a measure of the antiquark po- 
larization at x 2 in one proton, with the other proton providing a quark at X\ with known 
polarization |72|, f73|| . (This derivation also requires that the DIS structure function, gx(x), 
be approximately zero, a condition well met for x > 0.2.) 

A closely related reaction, the production of W ± , has been analyzed theoretically [[771 . 
The parity- violating nature of W ± production leads to some essential differences. First, only 
one of the two beams needs to be polarized. Second, W ± leptonic decay produces a charged 
lepton, which can be detected, and a neutrino, which cannot. Thus, the measurement 
necessarily integrates over some range of x\ and x 2 . Using the kinematics treated previously 
for DY production, one finds, for Ws produced at negative xp (opposite to the direction of 
the polarized beam (xi)), 

AT^-^-M, and AT^- — {x 2 ). (13) 

(X u 



2 Quarkonium production has also been considered as a means of determining gluon polarization. However, 
unresolved issues about the mechanism of quarkonium production (see Section 4) prevent its application as 
even a qualitative measure of AG, at least in the early stages of polarized hadronic collisions. 

3 Next-to- leading order QCD corrections to Equation |ll| have been evaluated by several authors [[74| [75| [T6| . 
They are generally small in most of the DY kinematic range. 
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This is similar to Equation but without the quark polarization factor multiplying the 
antiquark asymmetry. For Ws produced at positive xp, one finds, 

AT « — (x 2 ), and AT * ^(x 2 ). (14) 

In the DY process (or in Z° production, but not W ± production), it is also possible to 
measure a new structure function, called transversity, which is a correlation between quark 



momentum and its perpendicular spin component [7S]. The transversity is not measurable 
in inclusive DIS [f79|. It is measurable, in principle, in collisions of polarized protons whose 
spins are aligned perpendicular to the plane of dilepton detection. A nonzero transverse spin 
correlation in the DY process would clearly require both quark and antiquark transversities 
to be nonzero. A theoretically ideal, but currently impossible, experiment would be the 
measurement of the transverse spin correlation and hence the quark transversity of the 
proton via polarized pp collisions.^ 

Studies of both continuum DY and production in polarized p-p collisions are planned 



at RHIC 80|| . The scaling properties of the DY cross section, 

d 2 a/d yfrdxp oc 1/s, (15) 

combined with the kinematic relation, Equation |^, greatly favors the lowest beam energy 
consistent with the production of DY pairs in the "safe" region. Thus, detection of muon 
pairs with x\ = 0.25, x 2 = 0.4 (M = 5 GeV/c 2 ) at \/s = 50 GeV has 16 times higher cross 
section than detection of 20-GeV/c 2 pairs at yfs = 200 GeV. Because of its high threshold, 
W ± is feasible only at the highest proton energy of RHIC, y/s = 500 GeV. 

3.3 Charge Symmetry Violation in Parton Distributions 

Charge symmetry is believed to be well respected in strong interaction. Extensive exper- 
imental searches for charge symmetry violation (CSV) effects in various nuclear processes 
reveal an amount on the order of 1% [5T|. This is consistent with the expectation that CSV 



effects are caused by electromagnetic interaction and the small mass difference between the 
u and d quarks |82| . 



It has been generally assumed that the parton distributions in hadrons obey charge 
symmetry. This assumption enables one to relate the parton distributions in the proton and 
neutron in a simple fashion, u p (x) = d n (x),d p (x) = u n (x), etc. Indeed, charge symmetry is 
usually assumed in the analysis of DIS and DY experiments, which often use nuclear targets 
containing both protons and neutrons. Charge symmetry is also implicit in the derivation of 
many QCD sum rules, including the Gottfried sum rule, the Adler sum rule, and the Bjorken 
sum rule. 

The possibility that charge symmetry could be significantly violated at the parton level 
has been discussed recently by several authors |83], |84], ^5], |86|, |88|, jS9fl . Ma and collab- 



orators [RBI pi] pointed out that the violation of the GSR can be caused by CSV as well 



4 An estimate of the transversity asymmetry in the bag model ff9[ indicates that it peaks at large x, 
as does the helicity asymmetry. Polarized pp collisions at large x would thus yield a large experimental 
asymmetry. 



13 



as by flavor asymmetry of the nucleoli sea. They also showed that DY experiments, such 
as NA51 and E866, are subject to both flavor asymmetry and CSV effects. Using the bag 
model, Rodionov et al |^7|] showed that a significant CSV effect of ~ 5% could exist for the 



"minority valence quarks" [i.e. d p (x) and u n (x)] at large x (x > 0.4). A model study |[9| 



of CSV for sea quarks shows that the effect is very small, roughly a factor of 10 less than 
for valence quarks. Londergan & Thomas have recently reviewed the role of CSV for parton 
distributions 



Evidence for a surprisingly large CSV effect was recently reported by Boros et al |92| 
based on an analysis of F% structure functions determined from muon and neutrino DIS 
experiments. A large asymmetry, d n (x) ~ 1.2hu p (x) for 0.008 < x < 0.1, is apparently 
needed to bring the muon and neutrino DIS data into agreement. How would this finding, 
if confirmed by further studies, affect the E866 analysis of the flavor asymmetry? First, 
CSV alone could not account for the E866 data. In fact, an even larger amount of flavor 



asymmetry is required to compensate for the possible CSV effect |92| . Second, there has been 
no indication of CSV for x > 0.1. Thus, the large d/u asymmetry from E866 for x > 0.1 is 
not affected. 

3.4 Parton Distributions of Mesons and Hyperons 

Dilepton production using meson or hyperon beams offers a means of determining par- 
ton distributions of these unstable hadrons. Many important features of nucleon parton 
distributions, such as the flavor structure and the nature of the nonperturbative sea, find 
their counterparts in mesons and hyperons. Information about meson and hyperon parton 
structure could provide valuable new insight into nucleon parton distributions. Furthermore, 
certain aspects of the nucleon structure, such as the strange quark content of the nucleon, 
could be probed with kaon beams. 

Pion-induced DY cross sections have been measured in several high-statistics experi- 
ments |5B], [H], |§B|. These data form the basis for a global analysis |F7| to extract the 



pion structure functions. Although a large amount of DY data exists for tc~ beams, the 
corresponding data for 7r + beams are surprisingly meager. The tt + data are crucial for sep- 
arating the valence and sea-quark distributions in pions. The lack of high-statistics tt + DY 
data is responsible for our poor knowledge of the sea-quark distributions in pions. Future 
DY experiments using high-energy ir + beams (P n > 400 GeV/c) are required to study the 
sea of the pion. 

The advent of the Fermilab Main Injector (FMI) opens the possibility of performing DY 
measurements with intense pion and kaon beams. The relatively low beam momenta are 
suitable for studying parton distributions at large x. Londergan et al MB suggested that a 



comparison between n + and n DY cross sections on hydrogen and deuterium targets could 
test the charge symmetry of the valence quark distributions in the nucleon. 

Very little data exist for the kaon-induced DY process. The NA3 collaboration PS 



obtained several hundred DY events with a K~ beam. By comparing the K~ with the ti~ 
DY data, they found evidence that the u distribution in K~ is significantly softer than in 
the 7r~. Because this effect is the largest at large x±, the kaon beam at the FMI could be 
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used for further studies. 



No data exist for hyperon-induced dilepton production. The observation of a large d/u 
asymmetry in the proton has motivated Alberg et al |100|| to consider the sea-quark 
distributions in the S. The meson-cloud model implies a d/u asymmetry in the S + even 
larger than that of the proton. However, the opposite effect is expected from SU(3) sym- 
metry. Although relatively intense S + beams have been produced for recent experiments at 
Fermilab, this experiment appears to be very challenging because of large pion, kaon, and 
proton contaminations in the beam. 

3.5 Nuclear Dependence of the Drell-Yan Process 

The famous nuclear dependence of DIS, referred to as the European Muon Collaboration 

[1021 . 



(EMC) effect, was discovered in the early 1980s ||101 



Figure 11 shows the general 



features of the EMC effect, with data from the NMC (TU3 . 

In 1986, the Fermilab E772 collaboration proposed to measure the nuclear dependence 
of the DY process as a means of further elucidating the EMC effect and the issue of parton 
structure in nuclei in general. It is clear from the dominance of the term u{x-\)u{x2) in 
the kinematic region xf > 0.2, that the proton- induced DY process provides a view of the 
nucleon that is complementary to DIS. In the E772 configuration, typically more than 90% 
of the cross section arises from this term. In contrast, only about 15% of the cross section 
in DIS at x ~ 0.1 is due to antiquarks. 

The results from E772 ||104|| (Figures |TT] and [L2]) show that there is no enhancement of 
the antiquark distribution in nuclei. Even in the absence of a specific model calculation, 
the lack of an antiquark enhancement in nuclei seems at odds with the picture of nuclei as 
nucleons bound by the exchange of mesons. (Models of the EMC effect have been reviewed 
recently [ |102|| . We do not present a complete overview here.) 

In order to be more quantitative, it is necessary to invoke a model calculation. The 
pion excess model, based on the convolution framework of Sullivan |56| , is most specific in 
predicting the nuclear antiquark distribution. After many years of studying this issue (see 
e.g. [ |105| |), the conventional wisdom is that there are "excess" pions in nuclei, which are 
in part responsible for nuclear binding ||106| . In terms of the pion excess model, the E772 
results set stringent limits on the collectivity of the ttNN vertex inside nuclei. It appears 
that the nuclear pion field is not collective at all; there are no more pions surrounding the 
average nucleon in a heavy nucleus than there are in the weakly bound system, deuterium. 
This contradicts conventional wisdom and is also at odds with sophisticated new solutions 
of the nuclear many-body problem using realistic nuclear forces ||106| , |107| , |108| , |109|| . When 
these techniques are applied at central nuclear-matter density, they predict a substantially 
increased pion density, An n = 0.18 per nucleon [106 , |108|| .p| 



Using the Sullivan model, Brown et al | |110| propose a solution that appears to maintain 
much of the conventional understanding that nuclear binding results from meson exchange. 
They achieve agreement with the E772 antiquark ratio by postulating a decrease in the effec- 



ence [106 



Reference [108 contains a correction to the pion excess per nucleon calculated originally in Refer- 
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tive masses of hadrons inside nuclei resulting from a partial restoration of chiral symmetry. 
In view of the strong recent interest in chiral symmetry restoration in relativistic heavy-ion 
collisions, this is a fascinating explanation. However, it is far from being universally ac- 
cepted as the answer to the mystery of antiquarks in nuclei. More recently, Koltun ||111|| has 
advanced another suggestion that may reconcile the absence of an antiquark enhancement 
with conventional nuclear theory. 

3.6 Shadowing 

In the past decade, shadowing — the reduction of the cross section per nucleon for nuclear 
targets at small x — has been very well characterized experimentally in DIS [ 103 1 . The 



experimental signature appears clearly in Figure [II], where the DIS ratio falls below unity 
for x < 0.08. Theoretically, shadowing has been studied extensively in the past 10 years 
(see [102 1 for many seminal references). A recent study by Kopeliovich et al [112] presents 



an alternative view to the parton recombination picture of Mueller & Qiu ||113 



Shadowing is also expected in hadronic processes. However, the beam energies available 
at fixed-target facilities have limited investigations thus far. The reason is clear from the 
following numerical example. In order to have a large enough cross section, one needs 
x\ < 0.5. This limit, combined with the requirement M > 4 GeV/c 2 at ^ = 38.9 GeV (800 
GeV/c fixed target) in Equation |, yields x^ > 0.02. Similar arguments apply to fixed-target 
inclusive direct photon production. To date, the only experimental evidence for shadowing in 



hadronic reactions is the reduction in the nuclear dependence seen in Figures [11] and [12| from 
E772 |104|1 . The DY and DIS ratios for the same targets show similar behavior, although 



the x range is more limited for the DY data. 

It has been suggested that the reduction in the DY nuclear dependence ratio at small x<i 
is more appropriately ascribed to the effects of initial-state energy loss ||114| , 115|| . This is 



unlikely in view of Figures [12| and [13| and the discussion in Section 3.7. Figure [12] shows the 
nuclear dependence ratios for two different ranges of Q 2 (M 2 ). It is clear that the reduction 
at small x<i persists despite the signifcant change in the xp range for the different Q 2 cuts. 
An alternative view is obtained by fitting all four target ratios in terms of the common 
parameter a, defined by 



<?a = CiV x A c 



(16) 



A deviation of a from unity indicates a cross section that is not proportional to the number 
of target nucleons. Figure ^ shows a versus xp for two bins of X2, one in the shadowing 
region (■upper) and one outside it(lower). The lower plot exhibits no statistically significant 
dependence on xp or x\, again showing that the suppression is correlated with the target 
momentum fraction. 

Shadowing in the p + A DY process is largely due to antiquarks in the nucleus, unlike in 
DIS, where quarks dominate for most of the explored region. Although shadowing effects are 



expected for antiquarks and gluons, there is no known requirement that they be identical ||116 



to those for quarks. Nevertheless, direct comparison of DY and DIS data for Ca versus 
2 H (Figure [ll]) shows them to be identical within statistics. In contrast, in the region 
0.1 < x < 0.2, where DIS shows a modest enhancement in the ratio, no corresponding 
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increase occurs in the DY ratio. This is also evident in Figure n, which compares Fe versus 
2 H for DIS and the DY process. 

3.7 Multiple Scattering and Energy Loss 

The energy loss of fast partons moving through hadronic matter has received extensive 
theoretical attention [|114], |Tj| [119|, pO], [121 



Recently, Baier et al ||121|| have derived a 



connection between the longitudinal energy loss and the mean p t accumulated in multiple 
parton-nucleon scatterings that is directly applicable to experiment. Specifically, 



dE/dz = -a s N c pt, 



(17) 



where N c = 3. The theory of Baier et al yields the nonintuitive result that the total energy 
loss is proportional to square of the path length traversed. 

The nuclear dependence of the DY process offers a particularly clean way to test this 
relationship, since multiple scatterings are confined to the initial-state quark's traversal of 
nuclear matter. CERN experiment NA10 [ |122|| made the first measurements of the nuclear 
dependence of the DY process as a function of p t . The ratio a a/ &n was found to increase with 
increasing p t . E772 ||123|| extended the nuclear dependence measurements by determining 
A(pf), where A(pf) = (pf)(A) — (p 2 )( 2 H), for DY production for a series of nuclei. Table 
1 lists A(p 2 t ) from NA10 and E772. The E772 values (PL McGaughey, JM Moss, JC Peng, 
unpublished data) have changed slightly^ from those given in Reference ||123[| . Figure |T5 
shows a fit to the E772 data with a nuclear dependence of the form (pf)(A) oc A 1 / 3 . 

To convert A(pf) to an energy loss, one can approximate a nucleus by a uniform sphere 
of radius RqA 1 / 3 . This yields a mean path length of jRqA 1 / 3 for the initial-state partons. 
Then, using Equation [TJ|, one finds AE « 0.59 GeV for the Tungsten data from E772. In 
the lab frame, Sx i ~ AE/Eb eam = 7.4 x 10~ 4 , an exceedingly small nuclear effect versus x\ 
or xf- The dashed line in Figure [lower plot) is a linear fit to the data. This slope can 
be converted into an initial-state energy loss by using the shape of the accepted spectrum 
versus X\. This yields AE = 2.0 ±1.7 GeV, which is consistent with, but weaker than the 
limit determined from the p t broadening. 

The absence of a sizable quark energy loss is in qualitative accord with the finding of 
a weak nuclear dependence for the production of leading hadrons in DIS ||124| , |125|| . As 
has been noted before ||121|| , however, an interpretation of the nuclear dependence of the pt 
broadening in photon-induced dijet production [|126jl via Equation |17| leads to a very large 
energy- loss effect. We have no new insight into this dilemma. However, it is clear that further 



experimental tests of Equation [L7| are needed. Understanding how fast partons propagate 
through hadronic matter is crucial for understanding highly relativistic heavy-ion collisions 
and the possible production of the quark-gluon plasma. 

4. QUARKONIUM PRODUCTION 

4.1 Quarkonium Production in Hadronic Collisions 

6 The new evaluation is based on an analysis of da/dpj for 2 H 10 , not available previously. 
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Much of this section focuses on aspects of the nuclear dependence of quarkonium produc- 
tion. However, it is important to appreciate the impact of the experimental results from the 
Tevatron Collider on the theory of quarkonium production in hadronic collisions. Braaten 
et al have recently reviewed this subject. 

In 1995, the CDF collaboration published cross sections ||127|| for vertex-identified (i.e. 
not arising from 6-quark decay) J/ip and ip' production at very large p t — a particularly 
advantageous region for comparison with theory. The result was that the time-honored 
color-singlet production model underpredicted the cross section by factors as large as 30. 
Motivated by this spectacular failure, intense theoretical activity has focused on more general 
production mechanisms, including color-octet channels. Color-octet production has become 
the new paradigm in theoretical descriptions of quarkonium production. 

At fixed-target energies, however, the production of quarkonia at relatively low p t presents 
a challenge to the new models ||128| , |129| , |130|| . The effects of color-octet channels are reduced, 
and there are still problems to account for the observed small polarization. It is clear that 
the continuing theoretical debate about mechanisms of quarkonium production will influence, 
and be influenced by, the large nuclear dependence effects discussed below. 

4.2 Nuclear Dependence of Quarkonium Production 

In contrast to the DY process, large nuclear effects are found in the hadronic production 
of the J/ij) resonance [ |131| , |132| , |133| , |134|| . Whether induced by protons, antiprotons, or pions 



over a wide range of beam energies, there is a significant reduction in the cross section per 
nucleon for heavy nuclei. [] 



Using 800-GeV/c protons, E772 [ I57 | was able to resolve the J/ip and ip' resonances for 
the first time in a nuclear-dependence measurement. Surprisingly, the E772 collaboration 
found that the nuclear dependences for the two states were identical within statistics. Equally 
unexpected, they found a significant nuclear dependence in the yields of the T , T' , and T" 
resonances ||123|| , although less than that observed for the J/ip and ip' . Again, within worse 
statistics, the nuclear dependences of the T and the combined T' and T" were the same. 
Figure [To] summarizes the integrated nuclear-dependence data from E772. 



CERN experiment NA38 [|138j| has confirmed the observation of equal nuclear suppression 
of J/ip and ip' states mp+A collision at 200 GeV/c, employing a modified version of the NA10 
spectrometer. NA38 went further to record the first detection of charmonium production 
in relativistic heavy- ion collisions ||138|| . A qualitatively new effect was found: The ratio 



0>'/ <T J/V> decreased by nearly a factor of two in S + U collisions with respect to the ratio in 



p + A collisions. The NA38 collaboration and its successor NA50 |37] have developed the 
capability to measure high-mass dimuon events in coincidence with event charge multiplicity 
and calorimetric transverse energy. They have recently recorded interesting data on J/ip 
and ip' production in Pb + Pb collisions at 157 GeV/c/nucleon (See the review by Gerschel 



& Hiifner [139] in this volume. 



The high-statistics experiment NA3 ||134|| confirmed observations by earlier experiments 



7 There are also significant nuclear effects in J/ip production with real |135| and virtual photons |136[ . 
These are usually interpreted via the vector-dominance model in terms of J /ip> absorption on nuclcons. 
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of large changes in the nuclear dependence as functions of pt and xp. The NA3 nuclear 
dependence ratios versus xp and pt have similar shapes for beams of 7r~, 7r + , and protons, 



as shown in Figures |T7| and p~8| . 

Figures |T9| and |20| compare proton-induced J/ip production at 200 GeV/c and 800 GeV/c 
versus xp and x%- Figure clearly demonstrates the lack of scaling with x%. Such scaling 
would be expected for an effect whose origin was primarily due to a nuclear dependence of 
the target structure function — such as shadowing. In addition, even though the 800-GeV/c 
data are in the range where shadowing is expected, the observed nuclear dependence is much 
larger than expected for shadowing alone. 

The nuclear dependence shows a rise when plotted against p t for the J/ip (Figure 



and the T resonances. Table 2 lists A(p^) for J/ip and T production. The NA3 values are 



from Badier et al ||134j| . The E772 values (PL McGaughey, JM Moss, JC Peng, unpublished 
data) use a parameterization that reproduces da /dpi for the T data QTOj on 2 H, and the p t 



dependence of the ratio of cross sections R{ga/g2 H ). Unfortunately, the same analysis could 
not be performed for the E772's J/ip data owing to the very limited acceptance in p t . 

The values of A(p^) are similar for the J/ip and the T, albeit at different beam energies. 
The effect is three to four times larger for quarkonia than for DY production. A factor of 
two of this increase may be attributable to final-state multiple scattering, which is absent in 
the DY process. 

4.3 Nuclear Dependence of Open Charm Production 

The nuclear dependence of open heavy flavor production is much more difficult to measure 
than that of DY or quarkonium production. Four experiments employing vertex detectors 



to identify D meson decays have determined values of a (see Equation 16). As shown in 
Table 3, in contrast to J/ip and ip' production, there is no evidence of a nuclear dependence 
to open charm production at small xp. No measurements have yet been made of the nuclear 
dependence of open beauty production. 

4.4 Nuclear Effects in Quarkonium Production and the Quark-Gluon Plasma 

Much of the interest in the nuclear dependence of quarkonium production is connected 
to the search for production of the QGP in collisions of relativistic heavy ions (see Gerschel 
& Hufner ||139|| , this volume). In brief, the concept put forth by Matsui & Satz |144|1 is that 



states such as the J/ip and ip' would not survive if a QGP were formed, since their radii 
exceed the effective cc screening length in the QGP. Thus, one should observe a suppression 
of these states in central collisions of heavy ions compared with p + p or p + A collisions 
(where the QGP is presumably not formed). The large nuclear effects found in p+A collisions 
complicate this simple picture. In recent years, much effort has been directed at refining the 
phenomenology of the Matsui-Satz effect in order to identify truly new physics that may be 
manifested in hot-dense matter occupying an extended volume | |139| . 

To what extent are nuclear effects in quarkonium production understood? As skeptical 
experimentalists, we are tempted to respond, "not very well." We do not attempt a detailed 
review of the numerous models on the market but rather enumerate some of the basic ex- 
perimental elements that must be addressed by models of quarkonium production in nuclear 



19 



systems. 



1. Significant nuclear effects are found in hadronic production of all heavy quarkonium 
states. Charm is more suppressed than beauty, but both exhibit nuclear dependences 
that are much larger than those found for the DY process at comparable values of Q 2 . 
The few measurements of open charm production, based on vertex-identified decays of 
D mesons, exhibit no nuclear dependence. 

2. The nuclear suppression for J ftp and ip' is the same, within errors, as it is for T and 
the sum of T' and T" resonances. Thus, the effect does not depend on the final size of 
the hadron. It is important to appreciate that roughly half of the total production of 
the J/ip originates in decays from higher states, including the if)' 



3. The nuclear suppression becomes more pronounced at large xp or small x^. However, 
it does not scale with x 2 - The nuclear dependence versus xp is similar for production 
of the J /if) with 200-GeV/c and 800-GeV/c protons. 

4. A relative nuclear enhancement occurs at large p t . While this effect is also observed 
for DY production ||104 , |122|| , the rise with p t is much larger for the quarkonium states. 

is more 



The interpretation of the pt dependence in terms of parton energy loss ||121 



complicated than for DY production (Section 3.6) because both initial- and final-state 
effects can contribute. 



Although there is now a considerable body of data, much uncertainty remains in the 
phenomenological description of quarkonium production in nuclei. New nuclear-dependence 
measurements at negative values of Xp would rank high as a discriminant between current 
competing models. Such measurements are difficult at fixed-target facilities, but it is hoped 
that they will be carried out in the early years of RHIC. 

Hadronic production of quarkonium itself remains a field of intense theoretical interest 
Many of the desirable new experimental investigations relating to the importance of color- 
octet production, e.g. measurements of polarization, would also have an immediate impact 
on issues related to nuclear dependence. Such measurements would also have a significant 
impact on the understanding of quarkonium suppression as a signature of the creation of the 
QGP. 

5. FUTURE PROSPECTS 

The next decade will see a qualitative change in experimental capability to study dilepton 
production in hadronic processes. Two new hadron colliders, RHIC and LHC, will be in 
operation. Both will have the capability of accelerating and colliding a variety of beams — 
p + p, p + A, to A + A, where A can range to the heaviest nuclei. RHIC will operate in 
the range 50 < y/s < 500 GeV/nucleon. Starting in 2000, RHIC will also be able to collide 
polarized protons up to t/s = 500 GeV. Although specific experimental capabilities are still 
being developed at these facilities, areas of new opportunity related to the subjects reviewed 
here can be readily identified. 
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p + p collisions have never been carried out above the energy range of the CERN inter- 
secting storage rings. Thus, DY production at RHIC and LHC will produce new information 
for parton density distributions as well as for other issues in QCD that are accessible via 
this well-understood reaction. Extension of shadowing studies in hadronic processes requires 
the highest energy p + A collisions. The DY and direct photon production are the most 
readily interpreted processes for such investigations. The DY process in central collisions of 
ultra-relativistic heavy ions is the most theoretically tractable of the QCD processes in this 
exceedingly complex dynamical system. Although DY production is expected to take place 
prior to the putative transition to the QGP, there may be surprises. 

W ± production in p+p collisions, a first for hadron colliders, will give a direct indication of 



the flavor asymmetry of the proton, not requiring the assumption of charge symmetry ||145 
Additionally, new tests of the validity of charge symmetry might be performed via reactions 
such as p + d — > W ± ( fl92"fl , S Vigdor, unpublished). As discussed in Section 3.2, polarized 
proton collisions offer a wealth of new opportunities for the delineation of different aspects 
of polarized parton structure. The enormous range of Q 2 spanned from DY to and Z° 
production will give a unique window for observing the interplay between quark and gluon 
spins. 

RHIC and the heavy-ion program at the LHC are strongly focused on finding evidence of 
the QGP. A major element in this search will be quarkonium suppression. As we have seen, 
however, the origins of many strong nuclear effects observed in quarkonium production in 
p + A collisions remain mysterious. With the new capabilities of dilepton measurement in 
collider detectors, some of these mysteries may be investigated more easily than in closed- 
aperture fixed-target experiments. 

The new 120-GeV/c Fermilab Main Injector (FMI), scheduled to begin operation in 
2001, and the proposed 50-GeV/c Japanese Hadron Facility ||146|| present capabilities at 
the other end of the x scale. In contrast with the past 800-GeV/c fixed-target program at 
Fermilab, which could only operate in competition with the collider, the FMI program will 
run in parallel. Simple scaling rules show that large- x physics is often best carried out at 
comparatively low energies, where cross sections are largest. As we have seen, measurements 
of d/u and the nuclear dependence of the antiquark sea at high x require beams from the 
FMI P7f 



In conclusion, dilepton production in hadron-hadron interaction continues to be a very 
powerful tool for probing the parton structure of hadrons and nuclei, and for studying the 
dynamics of the strong interactions. We expect significant progress in these areas for the 
forseeable future, as the next generation of hadron colliders and high-intensity accelerators 
begin operation. 
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Table 1: Values of A(p 2 ) = (jpf) (A) — (pf) ( 2 H) for the DY process. For NA10, the mass range 
for muon pairs was, respectively, M > 4.35 GeV/c 2 at 140 GeV/c and M > 4.2 GeV/c 2 at 
286 GeV/c. The mass range of the T (8.5 GeV/c 2 < M < 11 GeV/c 2 ) was also excluded. 
For E772, the range was M > 4 GeV/c 2 , with the exclusion of 9 GeV/c 2 < M < 11 GeV/c 2 . 



Experiment 



NA10 [T22[ 



E772 11231 



Beam 



140 GeV/c ir- 



286 GeV/c tT 



800 GeV/c p 



A 



W 



W 



c 



Ca 



Fc 



W 



A(p 2 )(GeV 2 /c 2 



0.16 ± 0.03 ±0.03(syst.) 



0.15 ± 0.03 ±0.03(syst.) 



-0.011 ± 0.015 ±0.03(syst.) 



0.046 ± 0.012 ±0.03(syst.) 



0.053 ± 0.012 ±0.03(syst.) 



0.106 ± 0.016 ±0.03(syst.) 
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Table 2: Values of A(p 2 ) = (j$)(A) — (p 2 )( 2 H) for proton-induced production of the J/ip 
and T resonances. The E772 values were determined starting with a function, da jdp\ = 
C * (1 + {pt/Po) 2 )~ 6 , with p = 2.8 GeV/c, which fits the 2 H data for the T resonance. The 
NA3 value for the J/ip uses X H rather than 2 H. 



Experiment 


State 


Beam 


A 


A(p 2 )(GeV 2 /c 2 ) 


NA3 13J 




200 GeV/c 


W 


0.34 ±0.06 


E772 P3| 


T 


800 GeV/c 


c 


0.141 ±0.078 ±0.03(syst.) 








Ca 


0.200 ± 0.041 ±0.03(syst.) 








Fe 


0.315 ± 0.039 ±0.03(syst.) 








W 


0.410 ± 0.078 ±0.03(syst.) 
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Table 3: Nuclear dependence of open charm production via vertex-identified D meson decays. 



Experiment 


Beam 


(x F ) 


a 


WA82 pq 


340 GeV/c tt" 


0.24 


0.92 ±0.06 


E769 [141 


250 GeV/c tt" 


0.0 


1.0 ±0.05 ±0.02 


E789 |142| 


800 GeV/c p 


0.031 


1.02 ±0.03 ±0.02 


WA92 [[23(1 


350 GeV/c vr" 


0.18 


0.95 ±0.06 ±0.03 
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Figure 1: Combined dimuon mass spectrum from Fermilab E866 0: p+p and p+d collisions 
at 800 GeV/c. The shape of the continuum results from the mass dependence of the DY 
process folded with the acceptance of the spectrometer. 
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Figure 2: Feynman graphs for the electromagnetic processes (a) deeply inelastic lepton 
scattering and (b) the Drell-Yan process. 
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Figure 3: Proton-induced Drell-Yan production from experiments NA3 [§] (triangles) at 400 
GeV/c, E605 § (squares) at 800 GeV/c, and E772 (Q; PL McGaughey et al, unpublished 
data) (circles) at 800 GeV/c. The lines are absolute (no arbitrary normalization factor) next- 
to-leading order calculations for p + d collisions at 800 GeV/c using the CTEQ4M structure 



functions [|TT|. 
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Figure 4: Drell-Yan angular distribution from Fermilab E772 [0: p + Cu collisions at 800 
GeV/c. The dimuons cover the mass region 11 < M M + M - < 17 GeV/c 2 with —0.3 < xp < 0.8 
and ^ < 6 GeV/c. Mean values for p 4 , xp , and M are 1.4 GeV/c, 0.16, and 11.9 GeV/c 2 , 
respectively. The solid curve is a fit to the data with the form 1 + Xcos 2 8, where A is 
0.96 ± .04 ± .06. 




Figure 5: Feynman graphs for order a s corrections to the Drell-Yan process 
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Figure 7: The ratio a pd /2a pp of Drell-Yan cross sections vs. X2 from Fermilab E866 @. The 
curves are next-to-leading order calculations, weighted by acceptance, of the Drell-Yan cross 
section ratio using the CTEQ4M [|ll[] and MRS(R2) [|5(J parton distributions. In the lower 
CTEQ4M curve, d — u was set to to simulate a symmetric sea. 
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Figure 8: The ratio of d/u in the proton as a function of x extracted from the Fermilab 
E866 H cross section ratio. The curves are from various parton distributions. Also shown 
is the result from NA51 plotted as an open box. 
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Figure 9: Comparison of the E866 |] d-ti results at Q 2 = 54 GeV 2 /c 2 with the predictions 
of pion-cloud and chiral models as described in the text. The data from HERMES |52| are 
also shown. 
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Figure 10: F% — F£ as measured by NMC at Q 2 = 4 GeV 2 /c 2 compared with predictions 



based on the MRS(R2) |5(| and MRST []53| parameterizations. Also shown are the E866 
results f2|, evolved to Q 2 = 4 GeV 2 /c 2 , for the sea-quark contribution to — F%. For each 
prediction, the top (bottom) curve is the valence (sea) contribution and the middle curve is 
the sum of the two. 
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Figure 11: Comparison of the cross section per nucleon ratios (Jca/ a2 H f° r DY produc- 
tion |104| and DIS [ |T03fl . 



42 



16<Q <36 



36<Q 2 <196 GeV 2 /c 2 



1.1 
1.0 
0.9 

S1.1 

«1.0 
c 0.9 

>T 1.1 
10 
Q 09 

1.1 
1.0 
0.9 



1 1 1 1 

: c/ 2 h 


. C/ 2 H I « 




: 

— i — i — U — i — 




:Ca/ 2 H 

T -* i 


1 1 1 1 

. Ca/ 2 H 






— i — i — i — i — 

Fe/ 2 H 


1 1 1 1 

. Fe/ 2 H _ 




■I - 


— i 1 1 1 — 

:w/ 2 h 


■ i ■ i 

. W/ 2 H . 
IT " T . 




. T. 


;i s 1 1 - 



0.0 0.1 0.2 0.1 0.2 



Figure 12: Ratio of Drell-Yan cross sections per nucleoli for heavy nuclei to deuterium 
versus target momentum fraction from E772 [ 1 04 1 . The two columns show data for different 
bins of Q 2 = M 2 . 
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Figure 13: Nuclear dependence coefficient a for the Drell-Yan process ||104|| versus xp for 
(a) x 2 < 0.075, right scale, and (b) x 2 > 0.075, left scale. The thin solid lines show a = 1. 
The dashed line is a linear least-squares fit to the lower points. Also shown is the mean value 
of x\ for (b) 
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Figure 14: Comparison of the cross section per nucleon ratios ap e /a2 H for DY produc- 
tion |104| and DIS [ pTTH - 
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Figure 15: A(p^) = ($f){A) - (p^)( 2 H) versus A for the DY process from E772 ([[L23|],PL 
McGaughey, JM Moss, JC Peng, unpublished data). The solid curve corresponds to 
0.027((A/2) 1 / 3 - 1). 
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Figure 16: Ratios of heavy-nucleus to deuterium integrated yields per nucleon for 800 
GeV/c proton production of dimuons from the Drell- Yan process and from decays of the 
J/ip, ip', T(1S), and T(2S + 3S) states ||104|| . The short-dash and long-dash curves represent 
the approximate nuclear dependences for the bb and cc states, A 96 and A ' 92 , respectively. 
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Figure 17: Ratios of J/ip cross sections per nucleoli, platinum over hydrogen, versus xp for 
200 GeV/c 7r + , 7r~, and proton data from experiment NA3 ||134| . 
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Figure 18: Ratios of J/ip cross sections per nucleon, platinum over hydrogen, versus p t for 
200 GeV/c 7r + , 7T", and proton data from experiment NA3 ||134| . 
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Figure 19: a versus xp for proton-induced J/tp production at 200 GeV/c ||134|| and 800 
GeV/c |37|. 
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Figure 20: a versus x 2 for proton-induced J/ip production at 200 GeV/c [|134|| and 800 
GeV/c |37|. 



